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Abstract 

Direct  numerical  simulation  was  performed  on  two 
patient- specific  abdominal  aortic  aneurysms  (AAA) 
using  physiologically  realistic  flow  conditions.  The 
patients  have  AAA  with  diameters  of  approximately  5 
and  7  cm,  respectively.  The  blood  flow  hemodynamics 
are  shown  to  consist  of  large-scale  periodic  structures 
and  the  flow  transitions  from  laminar-to-turbulent  and 
back  to  laminar  during  a  pulse.  Transition  does  not  oc¬ 
cur  until  after  the  maximum  flow  rate  has  been  reached 
and  turbulence  persists  in  the  AAA  after  the  bulk  flow 
decelerates  to  a  laminar  condition.  For  both  cases,  a 
jet  of  blood  forms  at  the  AAA  throat  and  impinges 
against  the  arterial  wall  at  the  site  of  maximal  aortic 
diameter.  During  the  pulse  structures  shear  across  the 
AAA  causing  locally  high  wall  shear  stress  (WSS). 
It  is  shown  that  the  cases  have  significantly  different 
hemodynamics.  The  smaller  sized  AAA  has  faster 
moving  flow  structures,  higher  WSS,  higher  turbulent 
kinetic  energy,  and  the  blood  spirals  down  the  AAA 
toward  the  iliac  arteries.  In  the  larger  sized  AAA,  the 
blood  flow  channels  down  on  the  anterior  side  and  a 
large  recirculation  zone  forms  on  the  posterior  side. 
A  large  recirculation  zone  causes  low  WSS  which  is 
know  to  promote  atherosclerosis. 

1  Introduction 

The  vascular  condition  known  as  abdominal  aortic 
aneurysm  (AAA)  is  a  dilatation  of  at  least  1.5  times 
the  diameter  of  the  aortic  segment  between  the  re¬ 
nal  arteries  and  the  iliac  bifurcation.  The  aorta  is  the 
longest  artery  in  the  human  body  and  is  approximately 

2  —  2.4  cm  in  diameter.  The  AAA  dilatation  weakens 
the  arterial  wall  making  it  susceptible  to  rupture.  Most 
AAA  continue  to  enlarge  progressively  and  the  risk  of 
rupture  is  related  to  size.  To  prevent  rupture,  AAA 
are  repaired  when  they  meet  a  size  criteria  of  5.5  cm 
in  men  and  5  cm  in  women  and  since  the  mortality- 
associated  risk  of  rupture  is  less  than  the  risks  posed 
by  surgical  repair  for  patients  below  the  diameter  cri¬ 
teria,  most  patients  with  small  aneurysms  post-pone 
surgery;  however,  as  many  as  33%  of  ruptured  AAA 
have  diameters  smaller  than  5  cm.  Currently,  there  is 
no  truly  reliable  way  to  evaluate  susceptibility  of  par¬ 
ticular  AAA  to  rupture  and  diameter  is  one  of  many 


Figure  1:  Geometry  for  (a)  case  1  and  (b)  case  2.  Units  are 
in  cm. 


criteria  that  could  be  used  to  determine  when  to  oper¬ 
ate  on  an  AAA.  Discovery  of  these  evaluation  criteria 
would  vastly  improve  AAA  management  and  decrease 
mortality. 

AAA  formation  and  growth  occurs  due  to  a 
combination  of  biological  and  hemodynamic  factors 
(Lasheras  2010).  A  known  hemodynamic  factor  is  the 
frictional  drag  exerted  by  blood  flow  on  the  arterial 
wall,  referred  to  as  wall  shear  stresses  (WSS).  Stable 
unidirectional  WSS  is  associated  with  healthy,  intact 
vascular  endothelium,  whereas  low  flow  with  oscillat¬ 
ing  WSS  is  associated  with  atherosclerosis  (Tarbell  et 
al.  2014).  Evidence  also  suggests  that  high  WSS  with 
sharp  spacial  gradients  (Dolan  et  al.  2012)  can  dam¬ 
age  the  arterial  wall.  Low  oscillating  WSS  is  primar¬ 
ily  caused  by  zones  of  recirculating  blood  and  high 
WSS  by  blood  flow  impingement.  Since  the  abdom¬ 
inal  aorta  is  generally  a  uniform  cylindrical  is  shape 
in  most  cases  the  flow  would  remain  laminar.  How¬ 
ever  the  shape  of  an  AAA  could  cause  the  blood  flow 
to  transition-to-turbulence  increasing  the  oscillation  in 
the  flow  direction  and  causing  higher  WSS.  The  influ¬ 
ence  of  vascular  shape  on  turbulence  is  not  well  under¬ 
stood  and  the  goal  of  this  study  is  to  numerically  cal¬ 
culate  physiological  blood  flow  in  two  patient-specific 
AAA  to  determine  if  there  is  a  notable  different  in 
the  flow  structures  and  turbulence  intensity  between 
the  cases.  Previous  work  on  physiological  blood  flow 
in  patient- specific  arteries  includes  the  experimental 


study  by  Stamatopoulos  et  al.  (2011)  and  the  numer¬ 
ical  study  by  Les  et  al.  (2010),  Valen-Sendstad  et  al. 
(2011)  and  Arzania  &  Shaddenb  (2012). 

2  Computational  methodology 

The  AAA  geometry  used  in  this  study  is  generated 
from  CT  scans  using  the  commercial  medical  imaging 
software  Mimics.  The  geometry  of  the  two  patients, 
referred  to  as  case  1  and  case  2,  is  shown  in  figure  1 . 
Case  1  has  an  AAA  diameter  approaching  the  size  cri¬ 
terion  of  5  cm  and  case  2  has  a  AAA  diameter  that  has 
dilated  to  over  7  cm,  well  beyond  the  size  threshold. 
For  both  cases,  the  aorta  contorts  as  it  leads  into  the 
AAA.  To  simplify  the  geometry,  arteries  that  branch 
off  of  the  aorta,  such  as  the  renal  arteries,  are  excluded 
from  the  model  and  the  upstream  aorta  is  modeled  as  a 
constant  diameter  tube  that  transitions  into  the  patient- 
specific  aorta  shape.  The  front  and  back  sides  of  the 
AAA  are  referred  to  as  anterior  and  posterior,  respec¬ 
tively. 

The  period  of  a  blood  pulse  is  given  by  T  and 
the  phase  position  within  the  period  is  represented  as 
(j){t)  G  [0,T].  The  length  scale  is  the  upstream  aorta 
radius  R  and  velocity  components  are  normalized  by 
the  time  average  of  the  bulk  velocity  Tib,  where  an 
overbar  indicates  a  flow  quantity  has  been  averaged 
over  a  blood  pulse.  The  non-dimensional  parame¬ 
ters  that  govern  vascular  flow  are  the  bulk  Reynolds 
number  Re 5,  Womersley  number  Wo  and  the  max-to- 
mean  flow  rate  ratio  /3  defined  as: 


Re\)  =  2  u^Rjv^ 

(i) 

Wo  =  R\Jn/v, 

(2) 

—  ^b,max  / 

(3) 

Here  v  =  11/ p  is  the  kinematic  viscosity,  with  p  and  p 
being  the  dynamic  viscosity  and  density,  respectively, 
n  =  2tt/T  is  the  frequency  and  the  subscript  max  in¬ 
dicates  the  maximum  value  that  occurs  during  a  pulse. 
For  this  study,  using  flow  conditions  representative  of 
the  abdominal  aorta,  Re 5  =  390,  Wo  =  15.3  and 
/ 3  =  6.5.  Wo  describes  the  phase-lag  between  the 
pressure  gradient  and  the  flow  rate.  In  small  arter¬ 
ies,  where  Wo  «  1,  the  pressure  gradient  is  aligned 
with  the  flow  rate  and  it  can  be  treated  as  quasi-steady. 
As  Wo  is  increased,  the  flow  rate  begins  to  lag  be¬ 
hind  the  pressure  gradient  and  the  velocity  profile  no 
longer  matches  Poiseuille  flow.  Wo  >  10  is  consid¬ 
ered  highly  unsteady  with  inertia  forces  dominate  over 
viscous  forces  (Ku  1997).  The  parameters  used  in  this 
study  are  summarized  in  table  1 . 

An  incompressible  fluid  with  Newtonian  properties 
is  assumed.  The  governing  equations  for  the  conserva¬ 
tion  of  mass  and  momentum,  in  a  Cartesian  coordinate 
system  x  =  (x,  y,  z),  are  as  follows: 


Reb 

390 

Wo 

15.3 

P 

6.5 

P 

1050 

kg  m  3 

P 

0.0035 

kg  m-1  s' 

T 

1 

s 

R 

0.011 

m 

maximum  flow  rate 

147.2 

ml  s_1 

mean  flow  rate 

22.8 

ml  s_1 

Table  1 :  Summary  of  the  parameters  used  in  study. 


V  -  u  =  0  (4) 

du/dt  +  u  •  Vm  =  — Vp  +  z A72m  (5) 

where  the  velocity  field  is  given  by  u  =  (u,  v,  w),  p  = 
P/ p,  with  P  being  the  pressure. 

The  open  source  finite- volume  code  OpenFO AM- 
2.2.0  is  used  to  directly  solve  the  governing  equations. 
To  ensure  the  Courant-Friedrichs-Lewy  (CFL)  num¬ 
ber  <  1,  the  time  step  is  fixed  at  At  =  0.5  x  10-4  s 
during  a  pulse.  The  simulation  is  performed  on  48  pro¬ 
cessors.  The  initial  flow  field  is  set  to  zero  and  statis¬ 
tics  are  collected  after  10  pulses  have  been  simulated 
to  remove  the  effect  of  initial  conditions.  An  unstruc¬ 
tured  grid  composed  of  3078214  and  4462192  tetra¬ 
hedral  and  prism  cells  for  case  1  and  case  2,  respec¬ 
tively,  are  generated  using  the  commercial  grid  gener¬ 
ation  software  Pointwise.  The  wall  tangent  and  nor¬ 
mal  edge  length  of  the  first  prism  cell  adjacent  to  the 
wall  is  given  by  A st  =  0.005  cm  and  A sn  =  0.075 
cm,  respectively.  Using  the  maximum  WSS  at  the  in¬ 
let  boundary  to  scale  the  grid  spacing  in  viscous  units, 
A =  0.8  and  A =  11.6.  The  superscript  +  de¬ 
notes  scaling  with  the  viscous  length  v/ur,  where  ur 
is  the  friction  velocity.  From  the  wall  14  prism  cells 
are  extruded  at  an  expansion  rate  of  1.1.  The  tetra¬ 
hedral  cells  that  occupy  the  remaining  of  the  domain 
have  a  viscous  edge  length  of  As+  =  11.6. 

Systole  refers  to  when  blood  is  ejected  from  the 
heart  and  diastole  to  when  the  heart  muscle  relaxes. 
The  phase  variation  in  the  abdominal  aortic  flow  rate 
used  in  this  study  is  shown  in  figure  2  (a).  The  pulse 
profile  originated  from  Les  et.  al.  (2010)  as  discrete 
flow  rate  data  and  is  expressed  as  a  12  term  Fourier 
series.  A  fully  developed  unsteady  velocity  profile  is 
calculated  from  the  Fourier  coefficients  using  proce¬ 
dures  given  in  Womersley  (1955),  Hale  et  al.  (1955), 
McDonald  (1955)  and  Loudon  and  Tordesillas  (1998). 
This  profile  is  used  as  the  inlet  condition  upstream  in 
the  aorta.  Figure  2  (b)  shows  the  phase  variation  in 
WSS  at  the  inlet  boundary.  The  highest  WSS  during 
the  pulse  is  2.8  N  m-2  and  occurs  at  (f)/T  =  0.17. 
There  is  a  phase-lag  of  (j)/T  =  0.02  between  the  max¬ 
imum  WSS  and  maximum  flow  rate.  As  the  flow 
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Figure  2:  Phase  variation  at  the  inlet  boundary  in  (a)  Ub/ub 
and  (b)  WSS. 


decelerates  during  the  later  half  of  systole,  the  ad¬ 
verse  pressure  gradient  causes  the  flow  to  separate  at 
c j)/T  =  0.29  resulting  in  negative  WSS.  The  bulk  flow 
does  not  decelerate  to  zero  until  cj)/T  =  0.40  and  re¬ 
verses  in  direction  during  the  start  of  diastole.  The 
lowest  WSS  is  -1.1  N  m2  and  occurs  at  <p/T  =  0.40. 
During  the  later  half  of  diastole,  WSS  and  the  bulk 
flow  are  effectively  zero.  Typically  in  arteries  instan¬ 
taneous  WSS  ranges  from  1-7  N  m-2  (Malek  et  al. 
1999)  and  Dolan  et  al.  (2012)  classified  high  WSS  as 
>  3  N  m-2.  At  the  outlet  boundaries,  downstream  in 
the  iliac  arteries,  both  inflow  and  outflow  is  allowed 
depending  on  the  boundary  flux.  The  arterial  walls  are 
assumed  to  be  rigid  with  a  no-slip  condition. 

Phase  averaging  procedure 

Using  the  triple  decomposition  method  introduced 
by  Reynolds  &  Hussain  (1970)  a  periodic  flow  quan¬ 
tity  /(x,  t )  can  be  decomposed  into 

f(x,  t )  =  J(x)  +  fix,  <t>)  +  f\x,  t),  (6) 

where  /( x)  is  the  time  average,  f(x,  0)  is  the  peri¬ 
odic  component  and  f'(x,t )  is  the  turbulent  compo¬ 
nent.  Each  blood  pulse  is  segmented  into  M  —  25 
equidistant  phase  positions  and  at  each  phase  TV  =  50 
flow  fields  are  collected.  These  values  are  limited  by 
both  computational  and  storage  resources.  The  aver¬ 
age  of  /( x,  t )  locked  at  a  phase  is  calculated  as 

1  7V_1 

{f{x,ct>))  =  -Yjnx,t  +  iT),  (7) 
i= o 

with  (/  (je,  (/))  )  being  referred  to  as  the  phase  average 
and  is  the  sum  of  /( x)  and  /( jc,  0).  The  instantaneous 
turbulent  component  can  be  calculated  as 


Figure  3:  Case  1  isosurfaces  of  instantaneous  coherent  struc¬ 
tures  visualized  by  —20 ul/R2  using  the  A2  cri¬ 
terion:  phase  positions  (a)  (j)/T  =  0.28,  ( b ) 
0/T  =  0.32,  (c)  <j)/T  =  0.40,  (d)  (j)/T  =  0.56, 
(e)  <p/T  =  0.72  and  (/)  0/T  =  0.88. 


f(x,t)  =  f  (x,t)  -  (f  (x,4>))  (8) 

and  the  time  average  of  flow  quantities  is  calculated 
as  the  summation  of  the  phase- averaged  quantities, 
where 


m 


M- 1 

i=0 


(9) 


3  Results 

Figure  3  shows  isosurfaces  of  instantaneous  coher¬ 
ent  structures  for  case  1  visualized  by  —20 ul/R2  us¬ 
ing  the  A2  criterion  of  Jeong  &  Hussain  (1995).  Here 
A2  is  the  second  largest  eigenvalue  of  the  the  tensor 
S2  +  fi2,  where  S  and  are  the  symmetric  and  anti- 
symmetic  components  of  the  velocity  gradient,  respec¬ 
tively.  During  systole  flow  structures  form  at  the  AAA 
throat  and  by  mid  systole  impinge  against  the  arterial 
wall  at  the  site  of  maximal  aortic  diameter.  The  struc¬ 
tures  spiral  within  the  AAA,  shearing  against  the  ar¬ 
terial  wall  and  breaking  down  into  smaller  turbulent 
structures.  During  diastole  the  bulk  flow  through  the 
AAA  is  effectively  zero  but  the  turbulent  structures 
continue  to  spiral  making  approximately  1.5  revolu¬ 
tions  during  a  pulse.  Since  the  diameter  of  the  AAA  is 
larger  than  the  upstream  aorta,  the  local  Wo  is  higher, 


Figure  5:  Case  1  (tw,x)  along  the  intersection  of  the  arterial  wall  and  a  y-z  plane  at  x  =  12.5  cm:  phase  positions  (a) 
<p/T  =  0.24,  (6)  <j>/T  =  0.4,  (c)  <j>/T  =  0.56  and  (d)  </>/T  =  0.72. 


Figure  6:  Case  2  along  the  intersection  of  the  arterial  wall  and  a  y-z  plane  at  x  =14  cm:  phase  positions  (a)  4>/T  = 

0.24,  (b)  <j>/T  =  0.4,  (c)  tj>/T  =  0.56  and  (d)  4>/T  =  0.72. 


Figure  7:  Case  2  (tw,x)  along  the  intersection  of  the  arterial  wall  and  a  y-z  plane  at  x  =  21  cm:  phase  positions  (a)  <b/T  = 
0.24,  (b)  4>/T  =  0.4,  (c)  <j>/T  =  0.56  and  (d)  <j>/T  =  0.72. 


Figure  4:  Case  2  isosurfaces  of  instantaneous  coherent  struc¬ 
tures  visualized  by  —5 ul/R2  using  the  A2  cri¬ 
terion:  phase  positions  (a)  (j>/T  =  0.32,  ( b ) 
0/T  =  0.60  and  (c)  <f>/T  =  0.80. 


Figure  8:  Case  1  surface  contours  of  (a)  (rW;X),  ( b )  r'w^x 
and  (c)  TWjX,rms  at  phase  position  (j>/T  —  0.4. 


increasing  the  prevalence  of  unsteady  inertia  forces  in 
the  AAA.  At  the  start  of  diastole  due  to  the  reversal  in 
bulk  flow  direction  some  of  the  coherent  structures  in 
the  aorta  are  pulled  upstream.  By  the  end  of  diastole 
the  turbulent  structures  are  dissipating  and  the  process 
begins  again  with  the  next  pulse.  Figure  4  shows  iso¬ 
surfaces  of  instantaneous  coherent  structures  for  case 
2  visualized  by  —5 ul/R2  using  the  A2  criterion.  Sim¬ 
ilar  to  case  1 ,  flow  structures  form  at  the  AAA  throat 
but  it  is  not  until  early  diastole  that  the  jet  impinges 
against  the  arterial  wall  and  breaks  down  into  a  tur¬ 
bulent  flow.  During  diastole  the  turbulent  structures 
continue  to  flow  downstream,  shearing  along  the  ante¬ 
rior  side  of  the  AAA.  Since  the  AAA  is  larger,  relative 
to  case  1,  the  structures  are  slower  moving  and  it  is 
not  until  systole  of  the  next  pulse  that  the  structures 
are  ejected  from  the  AAA  into  the  iliac  arteries.  Simi¬ 
lar  to  case  1,  the  site  of  maximal  aortic  diameter  is  the 
location  the  structures  impinge. 

Figure  5  shows,  for  case  1,  the  phase-averaged  x- 
component  of  the  WSS,  along  the  intersection 

of  the  arterial  wall  and  a  y-z  plane  at  x  =  12.5  cm. 
The  anterior,  posterior,  left  and  right  sides  are  referred 
to  as  A,  P,  L  and  R,  respectively.  Phase  averaging  re¬ 
moves  the  instantaneous  incoherent  motion  leaving  the 
large-scale  coherent  structures.  At  4>/T  =  0.24  there 
is  a  sharp  spike  in  (rw,x)  on  the  posterior  side  below 
the  AAA  throat.  High  (rw,x)  only  occurs  briefly  at 
this  location  and  by  <\>/T  =  0.4  the  flow  structures 
have  spiraled  to  the  anterior  side  and  (rw^x)  has  de¬ 
creased  in  magnitude.  For  the  remainder  of  the  pulse 
(rw,x)  is  relatively  low  along  the  slice  and  fluctuates 
between  positive  and  negative  values.  Figure  6  shows 
for  case  2  (rw^x)  along  the  intersection  of  the  arterial 
wall  and  a  y-z  plane  at  x  =  14  cm.  A  sharp  spike  in 
(rWjX)  is  shown  to  occur  on  the  left  anterior  side  at 
4>/T  =  0.56.  Figure  7  shows  for  case  2  (rWjX)  along 
the  intersection  of  the  arterial  wall  and  a  y-z  plane  at 
x  =  21  cm.  This  slice  is  located  slightly  above  the  il¬ 
iac  arteries  and  here  (rW:X)  is  relatively  uniform  along 
the  arterial  wall.  This  indictates  that  the  large-scale 


(a) 


Figure  9:  Case  2  surface  contours  of  (a)  ( b )  r'wx 

and  (c)  Tw,x,rrns  at  phase  position  <j>/T  =  0.56. 


structures  that  form  at  the  AAA’s  throat  do  not  interact 
with  this  downstream  location.  Along  the  slice  (rWiX) 
changes  from  positive,  negative  or  «  zero  with  the  ex¬ 
pected  WSS  value  from  figure  2  (b). 

Surface  contours  of  WSS  are  shown  for  case  1  at 
(j)/T  =  0.4  and  case  2  at  (j)/T  =  0.56  in  figure  8  and 
9,  respectively.  It  can  be  seen  in  figure  8  (a)  and  figure 
9  (a)  large  structures  of  positive  and  negative  (r^) 
shear  across  the  AAA  with  the  magnitude  of  (rWjX) 
being  higher  for  case  1 .  Figure  8  (b)  and  figure  9  (b) 
show  the  instantaneous  incoherent  WSS  x-component 
r'w  x  .  For  both  cases  fluctuations  in  r'w  x  are  localized 
to  the  proximity  of  the  large  structures  and  follows  the 
structures  movement  as  they  shear  across  the  AAA. 
Figure  8  (c)  and  figure  9  (c)  show  the  phase-averaged 
quantity  rw^x^rrns,  where  the  subscript  rms  indicates 
the  operation  For  case  1,  rw^x^rrns  occupies 

a  larger  surface  area  of  the  AAA  and  is  at  a  higher 
intensity,  relative  to  case  2. 

The  phase-averaged  turbulent  kinetic  energy  is  de¬ 
fined  as  (TKE)  =  1/2 (uf  •  u').  Figure  10  shows  for 
case  1  (TKE)/ul  on  a  x-z  plane  at  y  =  5.3  cm.  High 
TKE  is  not  present  in  the  AAA  until  late  systole  and  it 


TKE )/ul 


Figure  10:  Case  1  contours  of  (TKE^/ul  on  a  x-z  plane  at 
y  —  5.3  cm:  phase  positions  (a)  <f>/T  —  0.36, 
\b)  (j)/T  =  0.44  and  (c)  <f>/T  =  0.52. 
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Figure  11:  Case  2  contours  of  (TKE)/ul  on  a  x-y  plane  at 
x  —  1.6  cm:  phase  positions  (a)  <j)/T  —  0.48, 
(6)  (j)/T  =  0.60  and  (c)  0/T  =  0.80. 


spirals  within  the  AAA  following  the  path  of  the  flow 
structures.  During  diastole  (TKE) /v%  dissipates  in  in¬ 
tensity.  Figure  11  shows  for  case  2  (TKE )/ul  on  a 
x-y  plane  at  x  =  1.6  cm.  High  TKE  is  not  present 
in  the  AAA  until  early  diastole  and  similar  to  case  1 
the  TKE  follows  the  path  of  the  flow  structures  and 
dissipates  during  late  diastole.  Case  2  has  a  zone  of 
(TKE) /v%  on  its  posterior  side  that  persists  at  a  low 
intensity  throughout  diastole.  For  both  cases  there  is 
a  phase-lag  between  the  highest  flow  rate  and  the  for¬ 
mation  of  TKE,  but  this  phase-lag  is  more  profound 
for  case  2. 

The  turbulent  component  in  the  flow  is  further  in¬ 
vestigated  by  viewing  the  velocity  signal  at  discrete 
locations  in  the  AAA.  For  case  1  VI  =  (12,4,  2) cm 
and  V2  =  (12,  5,  2) cm  are  probes  located  within  the 
AAA’s  bulge.  Figure  12  shows  the  temporal  variation 
in  u/ub  at  these  probe  locations  during  the  last  three 
simulated  pulses.  At  both  probe  locations,  there  is  a 
distinct  coherent  component  to  the  velocity  signal  that 
repeats  every  period  and  an  incoherent  or  turbulent 
component.  The  turbulent  component  at  V2  qualita¬ 
tively  appears  to  be  at  a  higher  intensity,  relative  lo  VI. 
Figure  13  shows  for  case  2  u/ub  at  V3  =  (15, 8,  2)cm 
and  V4:  =  (18,  5,  2) cm  during  the  last  three  simulated 
pulses.  V3  is  located  on  the  anterior  side  in  the  path  of 
the  flow  structures  and  V4,  is  located  on  the  posterior 


Figure  12:  Case  1  temporal  variation  in  u/ub  at  probe  lo¬ 
cations  (a)  VI  and  (b)  V2  during  the  last  three 
simulated  pulses. 
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Figure  13:  Case  2  temporal  variation  in  u/ub  at  probe  lo¬ 
cations  (a)  V3  and  (b)  V4  during  the  last  three 
simulated  pulses. 


side.  It  can  be  seen  qualitatively  at  V3  the  turbulent 
component  is  at  a  lower  intensity  and  frequency,  rel¬ 
ative  to  the  case  1  probes.  At  V 4  there  is  no  distinct 
coherent  motion  to  the  velocity  signal  and  qualitatively 
it  appears  to  be  dominated  by  low-frequency  turbulent 
fluctuations.  In  figure  14  for  case  2  at  V3  the  velocity 
is  decomposed  into  its  instantaneous  turbulent  com¬ 
ponent  u' /ub  and  plotted  over  the  last  10  simulated 
pulses.  By  viewing  v! /ub  it  can  not  be  determined 
where  in  the  period  of  the  pulse  the  highest  turbulent 
intensity  would  occur.  In  figure  15  phase  averaging 
has  been  applied  to  the  instantaneous  turbulent  com¬ 
ponent  and  (TKE )/ul  has  been  calculated  at  V3.  It 
can  be  seen  that  the  highest  (TKE) /v%  does  not  occur 
until  late  diastole  and  is  «  zero  during  systole. 

The  time-averaged  streamlines  are  shown  for  case 
1  in  figure  16  (a)  and  case  2  in  figure  16  (b).  Time 
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Figure  14:  Case  2  temporal  variation  in  u' /ub  at  probe  loca¬ 
tion  VZ  during  the  last  10  simulated  pulses. 
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Figure  15:  Case  2  phase  variation  in  (TKE) /ul  at  probe  lo¬ 
cation  V3. 


averaging  removes  the  time-dependent  flow  structures 
and  gives  an  overall  description  of  the  flow  field  and 
the  hemodynamic  forces  acting  on  the  arterial  wall. 
For  case  1  the  time-averaged  flow  field  consists  of  a 
jet  that  impinges  against  the  arterial  wall  and  spirals 
down  the  length  of  the  AAA.  There  is  no  zone  of  sep¬ 
arated  recirculating  flow  for  case  1 .  The  time-averaged 
streamlines  for  case  2  show  that  the  flow  also  impinges 
against  the  arterial  wall,  however  instead  of  spiraling 
like  case  1,  the  blood  flows  downstream  on  the  ante¬ 
rior  side  and  a  large  zone  of  recirculating  blood  forms 
on  the  posterior  side.  For  both  cases  the  AAA  bulges 
at  the  location  the  time-averaged  streamlines  impinge. 

The  time-averaged  TKE  is  shown  for  case  1  on  a 
x-z  plane  at  y  =  5.3  cm  in  figure  15  (a)  and  case  2 
on  a  x-y  plane  at  x  =  1.6  cm  in  figure  15  (b).  It  can 


Figure  16:  Time-averaged  streamlines  for  (a)  case  1  and  ( b ) 
case  2. 


Figure  17:  Contours  of  TKE/ul  for  (a)  case  1  on  a  x-z  plane 
at  y/R  =  4.8  and  ( b )  case  2  on  a  x-y  plane  at 
x/R=  1.4. 


be  seen  that  for  case  1,  high  TKE/u \  occupies  a  ma¬ 
jority  of  the  AAA’s  volume  and  the  highest  TKE/ul 
occurs  along  the  centre  axis  of  the  AAA.  For  case  2  a 
zone  of  high  TKE/ul  occurs  on  the  anterior  side  in  the 
AAA’s  bulge  and  occupies  a  relatively  small  amount 
of  the  total  AAA  volume.  A  smaller  zone  of  lower  in¬ 
tensity  TKE/ul  occurs  on  the  posterior  side.  Overall 
TKE/ul  f°r  case  1  h  is  significantly  higher  and  occu¬ 
pies  at  larger  amount  of  the  total  AAA  volume,  relative 
to  case  2. 

4  Conclusions 

Direct  numerical  simulation  has  been  performed 
on  two  patient- specific  AAA.  The  hemodynamics  vary 
considerably  between  the  cases  and  the  flow  is  found 
to  be  comprised  of  complex  periodic  and  turbulent 
structures.  To  analyze  the  phase- variation  during  a 
pulse  the  velocity  was  decomposed  into  it’s  instanta¬ 
neous  turbulent  and  phase-averaged  component.  For 
both  cases  the  blood  flow  transitions  from  laminar-to- 
turbulent  and  back  to  laminar  during  a  pulse.  There 
is  considerable  phase-lag  between  the  maximum  flow 
rate  and  the  formation  of  turbulent  flow  structures  in 
the  AAA.  During  systole,  for  both  cases,  a  jet  forms 
at  the  AAA’s  throat  and  impinges  against  the  arterial 
wall  with  the  AAA  appearing  to  bulge  at  the  impinge¬ 
ment  location.  During  diastole,  when  the  bulk  flow 
through  the  AAA  is  «  zero,  the  velocity  in  the  AAA 
is  not  zero  and  the  flow  structures  continue  to  shear 
across  the  arterial  wall.  High  WSS  and  TKE  occurs 
in  the  proximity  of  the  flow  structures.  The  following 
differences  exist  between  the  two  cases: 

•  For  case  1  the  flow  transitions-to-turbulence  dur¬ 
ing  late  systole  and  for  case  2  during  early  di¬ 
astole.  For  both  cases  the  turbulence  has  mostly 
dissipated  by  early  systole  of  the  next  pulse. 

•  Overall  the  flow  velocity  is  higher  in  case  l’s 
AAA.  This  results  in  faster  moving  flow  struc¬ 
tures,  higher  WSS  and  TKE.  High  TKE  occupies 
a  larger  amount  of  the  AAA  volume  for  case  1 . 


•  For  case  1  the  blood  flow  spirals  within  the  AAA 
and  for  the  time-averaged  flow  there  is  no  zone 
of  separated  recirculating  blood.  For  case  2  the 
blood  flows  downstream  on  the  anterior  side  and 
a  large  zone  of  recirculating  blood  forms  on  the 
posterior  side. 

Within  the  scope  of  the  two  AAA  studied  here 
the  moderately  sized  AAA  has  higher  turbulence  and 
WSS.  However  the  larger- sized  AAA  has  a  large  recir¬ 
culation  zone  which  is  not  present  in  the  moderately- 
sized  AAA.  Physiologically  both  high  turbulence  and 
a  large  recirculation  zone  could  damage  the  arterial 
wall. 
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